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ABSTRACT 

A mathematical model is developed to determine the heat 
transfer rate of a cylindrical condenser section of a rotating 
heat pipe. This model is coupled to an existing code and an 
analysis is accomplished on both a smooth and axially finned 
condenser. The results of this analysis are compared to those 
of a similar analysis performed on a tapered condenser heat 
pipe using identical geometric and operating parameters. 

Results of the comparison indicate cylindrical condensers 
have a lower heat transfer rate than equivalent tapered con- 
densers. This reduction in heat transfer rate is due to a 
greater condensate film thickness and is most significant in 
a smooth condenser. 

Axially finned condensers with triangular and rectangular 
fin profiles are also compared. The rectangular fins are 
assumed to have adiabatic tips. Results indicate the heat 
transfer rates for these two profiles vary by only 0.40 per 
cent for both tapered and cylindrical condenser. 
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I. INTRODUCTION 



A. THE ROTATING HEAT PIPE 

The rotating, wickless heat pipe is a closed container 
designed to transfer large amounts of heat from rotating 
machinery components. Essentially, it consists of three 
main components: a cylindrical evaporator section, a con- 

denser section which may be either tapered or cylindrical in 
shape, and a fixed amount of working fluid. A typical tapered 
rotating heat pipe is shown in Figure 1. 

When the heat pipe is rotated about its longitudinal axis 
at a speed above a certain critical value, the working fluid 
forms an annulus in the evaporator section. Note in Figure 1 
that the diameter of the evaporator is larger than the con- 
denser. This larger diameter provides a greater liquid reser- 
voir. As heat is added to the evaporator, the fluid in the 
evaporator will vaporize. The vapor will flow axially towards 
the condenser as a result of a slight pressure difference, 
transporting the latent heat of vaporization with it. In the 
condenser end, external cooling of the condenser causes the 
vapor to condense. In the case of a tapered heat pipe, the 
centrifugal force due to the rotation of the pipe has a com- 
ponent acting along the condenser wall which accelerates the 
liquid condensate back to the evaporator to complete the cycle. 
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Figure 1. Schematic Drawing of a Tapered 
Condenser Rotating Heat Pipe 
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A cylindrical heat pipe, on the other hand relies on a hydro- 
static pressure gradient to drive the liquid condensate back 
to the evaporator. 

B . BACKGROUND 

The first theoretical investigation into the performance 
of a tapered rotating heat pipe at the Naval Postgraduate 
School was accomplished by Ballback [Ref. 1] in 1969. He 
examined the limits in heat transfer controlled primarily by 
fluid dynamic considerations. In particular, he considered 
the following four limits on heat pipe performance: a) the 

boiling limit, b) the entrainment limit, c) the sonic limit 
and d) the condensing limit. Tantrakul [Ref. 2] calculated 
these limits for a specific heat pipe. He found the condens- 
ing limit was the controlling limitation. In fact, the cal- 
culated heat transfer rate, based on the condensing limit was 
l/10th the heat transfer rate for the next lowest limit, the 
entrainment limit. 

In order to overcome this condensing limitation and thus 
increase the heat transfer rate of the rotating heat pipe, 
the concept of an internally finned tapered rotating heat 
pipe was considered by Schafer [Ref. 3]. Schafer developed 
an analytical model for this tapered heat pipe with a trian- 
gular fin profile as shown in Figure 2. He assumed one dimen 
sional heat conduction through the wall and fin. Corley [Ref 
4] developed a two-dimensional heat conduction model using a 
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Figure 2. Axially Finned Condenser Geometry Showing 
Fins, Troughs, and Lines of Symmetry 
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Finite Element Method formulation for the same geometry. To 
overcome the problem of few nodal points along the fin surface, 
Corley assumed a parabolic temperature profile along the fin 
surface. Tantrakul [Ref. 2] modified Corley's computer code by 
increasing the number of finite elements from two to three in 
order to minimize the error at the apex of the fin. Purnomo 
[Ref. 5] developed a two-dimensional Finite Element .Method 
solution to the steady state heat conduction problem using a 
linear triangular finite element. Davis [Ref. 6] modified 
Purnomo ' s code to make it compatible with COPES/CONMIN [Ref. 8], 
an optimization program. Davis, in his modification, found a 
coding error in Purnomo's [Ref. 5] code, that once corrected, 
permitted Purnomo's corrected code to converge to Schafer's 
[Ref. 3] results. 

Purnomo's [Ref. 5] code is limited in that it is restricted 
to one particular condenser geometric configuration, namely: 
an axially finned tapered condenser heat pipe with a triangular 
fin profile. In that tapered finned condenser are difficult 
to manufacture, it is doubtful that any widespread practical 
application of this geometric configuration will result. Cylin- 
drical condensers, on the other hand, can be manufactured with 
much less difficulty and might find practical application. 

This being the case, a more practical and beneficial code 
would be one that could analyze cylindrical condenser rotating 
heat pipes, both finned and smooth. In actuality, the most 
beneficial code would be one that could analyze the following 
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four geometric configurations: 1) tapered-internally finned, 

2) tapered- smooth , 3) cylindrical-internally finned, and 4) 
cylindrical-smooth. The theoretical heat transfer performance 
of the four geometries could then be compared to determine the 
advantages and disadvantages of each design. An additional 
advantage would be gained if different fin profiles, i.e., 
triangular vice rectangular, could also be analyzed and 
compared . 

C. THESIS OBJECTIVES 

The objectives of this thesis are: 

1) Develop analytical models for both cylindrical - smooth 
and cylindrical -axially finned condensers. 

2) Develop solution techniques to these analytical models 
that will account for temperature variations along the 
axial length of the condenser. 

3) Modify Purnomo ' s [Ref. 5] code to provide a solution to 
the two-dimensional steady state conduction heat transfer 
problem for the following four geometric configurations: 
a) tapered-smooth , b) tapered- finned , c) cylindrical- 
smooth, and d) cylindrical-finned. 

4) Modify Purnomo ' s [Ref. 5] code to provide the additional 
capability of analyzing a rectangular fin profile with 
an adiabatic tip. 

5) Obtain and compare results of the four geometric con- 
figurations given above for various operating conditions. 



20 



II. THEORETICAL ANALYSIS FOR A CYLINDRICAL HEAT PIPE 



A. INTRODUCTION 

In a cylindrical condenser heat pipe, the radius of the 
condenser is constant along the axial length of the condenser. 
The flow of the condensate in the absence of vapor- liquid 
interfacial shear, is dependent upon the variation in hydro- 
static pressure with changes in film thickness along the sur- 
face of the heat pipe. Leppert and Nimmo [Refs. 8 and 9] 
investigated the phenomenon of film condensation on a flat 
horizontal plate. This situation is similar to film conden- 
sation on the inside surface of a rotating cylindrical con- 
denser. In the case of a cylindrical condenser, the body 
force, rather than being the force of gravity, is now the 
centrifugal force caused by the rotation of the heat pipe. 
Weigenseil [Ref. 10] and Tantrakul [Ref. 2] compared experi- 
mental results for a cylindrical condenser rotating heat pipe 
with the theoretical results of Leppert and Nimmo [Refs. 8 
and 9] and found good agreement. The Leppert and Nimmo solu- 
tion was limited in that it was based on a constant surface 
temperature along the length of the plate. A rotating heat 
pipe, in actuality, has a temperature variation along the 
axial length of the condenser which in some cases, may be 
significant. This being the case, it was necessary to develop 
a mathematical model which would consider the axial temperature 
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variation in the solution of the heat transfer analysis. In 
the mathematical development that follows, a cylindrical 
smooth (unfinned) condenser will first be considered in that 
it is the simplest case. The model will then be extended to 
include a cylindrical axially finned condenser. 

B. THEORY FOR A CYLINDRICAL SMOOTH CONDENSER 
1 . Assumptions 

In developing the theoretical analysis, the following 
assumptions are made: 

a) Film condensation, not dropwise condensation occurs in 
the condenser. 

b) The condensate film undergoes laminar flow. 

c) Momentum changes through the condensate are small. 

Thus, there is essentially a static balance of forces. 

d) The vapor exerts no drag in the condensate; there is 
no interfacial shear. 

e) The temperature distribution within the film is linear. 

f) The vapor space is essentially at one pressure, P . 

g) The density of the fluid is much greater than the den- 
sity of the vapor. Thus, the density of the vapor can 
be neglected. 

h) The centrifugal force is much greater than the force 
of gravity and, thus, gravity may be neglected. 

i) Velocity gradients in the circumferential direction 
relative to the pipe wall are negligible. 
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j) The condensate film thickness is much less than the 
radius of curvature of the condenser wall. 

k) The rotating heat pipe is operating at steady state 
conditions . 

2 . Condensate Momentum Equation (X-Direct ion) 

By applying the above assumptions and the coordinate 
system shown in Figure 3, an analysis similar to Nusselt's 
original film condensation film theory may be used [Ref. 11]. 
Based on assumption c, a static force balance may be taken 
on an infinitesimal fluid element in the x-direction as shown 
in Figure 3. This force balance results in the following 
equation: 

V 

ZF x = 0 ; T57 - H - 0 ^ 2- 1 ) 

where t = shear stress (lbf/ft 2 ) 
p = pressure (lbf/ft 2 ) 

x = co-ordinate measuring distance along surface (ft) . 
y = co-ordinate measuring distance normal to surface (ft). 

3 . Condensate Momentum Equation (Y-Direction) 

In a similar manner, using Figure 3, a force balance 
in the y-direction yields: 

ZF = ° ; + p w 2 r = 0 (eqn 2.2) 

y 3y 1 

where = density of the fluid (lbm/ft 3 ) 
go = angular velocity (rad/hr) 
r = radius (ft) 
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Figure 3. 



Cross Section of Infinitesimal Fluid Element on 
Cylindrical Condenser Internal Surface 
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4 . Fluid Velocity 

Integrating equation (2.2) between the limits y and 
6* for y and corresponding limits of P and P y for pressure 
results in the following equation: 

P = P y + P£ 0 ) z r (6*-y) (eqn 2.3) 



where P v = the pressure of the vapor (lbf/ft 2 ) , and 
6* = film thickness (ft) 

Differentiating equation (2.3) with respect to x yields the 
following expression for dP/dx: 



dP _ dP v + p£(D 2 r dS* 
dx dx dx 



( 

(eqn 2.4) 



Applying assumption (f) , (P y is constant, therefore, dP v /dx=0) 
and substituting equation (2.4) into equation (2.1) yields: 



3t 

3y 



p f u 2 r 



d6* 

dx 



(eqn 2.5) 



Integrating equation (2.5) with the corresponding limits of 
integration y to 6* and x to 0 results in the following ex- 
pression for shear stress: 



T 



p £ oo 2 r 



d6* 

dx 



[y - 6 * ] 



(eqn 2.6) 
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But , 



T 




(eqn 2.7) 



where y = fluid dynamic viscosity (lbm/ft-hr) 
u = condensate velocity (ft/hr) 

Substituting equation (2.7) into equation (2.6) and integrating 
with the corresponding limits of integration 0 to y and 0 to u 
yields : 



The average velocity of the condensate may be found in the 
following manner: 



p f u) 2 r d6 * 




(eqn 2.8) 



u 



dx 




(eqn 2.9) 



or 



u 




(eqn 2.10) 



5 . Continuity Equation 

The continuity equation for mass flow requires that: 



m 




(eqn 2.11) 



where m = condensate mass flow rate (lbm/hr) 

A = cross sectional area of the fluid (ft 2 ) 
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This can also be written as 




p^u 2iTrdy 



(eqn 2.12) 



Substituting equation (2. ID) into equation (2.12) and inte- 
grating yields : 



2 

2irpf io 2 r 2 d6* 
\i dx 




(eqn 2.13) 



Differentiating this equation with respect to x yields: 



dm _ _ 27rpf oj 2 r 2 d r dg_* g * 3 
dx U dx dx 3 J 



(eqn 2.14) 

v 



6 . Energy Equation 

Having applied assumption (e) , if the film surface 
temperature is at the saturation temperature (T sat ) of the 
vapor and if the wall of the axial increment is at a given 
constant temperature (T w ) , then the heat transfer by conduc- 
tion of a fluid element of surface area dA is: 



dq = C£sat ~ T w ) dA 
g " 5 



(eqn 2.15) 



where dq 
dA 

kj- 



differential heat transfer rate (Btu/hr) 
2rrrdx(f t 2 ) 

thermal conductivity of the condensate film 
(Btu/hr- ft-F) 
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T sat = saturation temperature (degrees F) 

T = inside condenser wall temperature (degrees F) 
w 

Considering the change of phase and defining K f as 
the average enthalpy change of the vapor in condensing to a 
liquid and subcooling to the average liquid temperature of 
the film, then dq is also defined by: 

dq = h £o . dm (eqn 2.16) 



where h £cT = laten heat of vaporization (Btu/lbm) 

Cp = specific heat (Btu/lbm R) 

AT = (T . - T ) 

v sat w J 

hr- = hr + 0.35 • c ♦ AT 
fg fg P 

Rearranging equation (2.16) and substituting this equation into 
equation (2.15) yields: 



dm _ k-f(T sa t- ~ T w ) 2TTr 
dx 7 s 



(eqn 2.17) 



Finally coupling the energy and continuity equations result in 
the following differential equation: 



6 * d_ r d6* 
dx dx 



A 3 

6 ] 



5kf(T.sat. ~ T w ) P 
p f 2 m 2 r E £g 



(eqn 2.18) 



Equation (2.18) can be solved using the Finite Element 
Method to provide the film thickness profile along the axial 
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length of a cylindrical condenser. Appendix A provides a 
detailed description of this solution. Once the film pro- 
file is known, a steady state two-dimensional heat conduc- 
tion analysis can be performed. 

7 . Determination of Heat Transfer Rate 

Assume that the cylindrical condenser section of the 
rotating heat pipe is divided axially into a number of incre- 
ments. Then for any axial increment of a cylindrical conden- 
ser, the differential heat flux can be determined by the 
following expression. 



dq 11 = 



(Tsat " Too) 

— 55 

5 thick ^ 1 

FT T E . 

f w ext 



(eqn 2.19) 



where Too = ambient temperature (degrees F) 

thick = thickness of the condenser wall (ft) 
k w = thermal conductivity of the wall material 
(Btu/hr- ft-F) 

h ext = external heat transfer coefficient 
(Btu/hr-ft 2 -F) 

Note the three terms in the denominator are the thermal re- 
sistances of the film, wall and external convection respectively. 

The differential heat transfer rate for any increment 
can be found by the following relationships: 

dq = dq 11 • Zrrdx (eqn 2.20) 
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or 



dq 



2ttt (T sa t - 

95 




thick 




Tw) dx 

+ TT” — 

ext 



(eqn 2.21) 



Equation (2.21) represents the total heat transfer rate for 
an incremental section of width dx. To find the total heat 
transfer rate for the entire cylindrical condenser, the in- 
cremental heat rates must be summed over the entire length 
of the condenser. Therefore: 



Q 



total 



NDIV 
Z dq 
i = l 



(eqn 2.22) 



where NDIV = total number of axial increments. 

C. THEORY FOR A CYLINDRICAL AXIALLY FINNED CONDENSER 
1 . Assumptions 

Referring to Figure 4, it is obvious that the analysis 
of a cylindrical internally finned condenser is more compli- 
cated due to the mass flow from the fins into the trough 
region between the fins. For this reason, in addition to 
the simplifying assumptions made for the smooth condenser 
which are listed in the previous section, the following 
assumptions must also be made: 

a) Referring to Figure 5, the mass flow along the fin sur- 
face does not flow axially in the x-direction, but only 
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along the surface of the fin in the z-direction into the 
trough. Thus, mass flow in the axial direction is only- 
permitted in the trough region between the fins. This 
is a reasonable assumption in that the film thickness 
along the fin surface is very small in relation to the 
film thickness in the trough. This being the case, the 
hydrostatic force in the x-direction on a fin fluid ele- 
ment will be much less than the centrifugal force com- 
ponent in the z-direction on that same fluid element 
forcing that fluid element into the trough. 

b) Just as in the axial direction, there is no pressure 
change along the surface of the fin in the z-direction. 

c) It will be assumed that the temperature along the con- 
vective surface of the fin is at a constant value (T aycr ) . 
This average fin surface temperature is the arithmetic 
average of the fin tip temperature and the fin base sur- 
face temperature where the fin intersects with the wall 
of the condenser. This is a valid assumption if the fin 
section is divided into a sufficient number of finite 
elements. Purnomo’s [Ref. 5] results indicate a less 
than one degree variation, even for very large fin half 
angles. This variation in temperature will have an in- 
significant effect on film thickness along the surface 
of the fin and can be neglected by using an average 
value . 
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Figure 4. Trough Section of Axially Finned Condenser 
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Figure 5. Axially Finned Condenser Symmetric Fin Section 
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2. Mass Flow in the X-Direction 



As a result of assumption (a) , the resulting momentum 
equation in both the x and y directions as well as the equa- 
tions for velocity and mean velocity are identical to those 
developed in section B for the smooth condenser and will not 
be redeveloped here. Looking now at mass flow in the x-direc- 
tion which is limited to flow only in the trough, the mass 
flow rate is given by the following expression: 



™to tal 



P f 2 u 2 r d6* 
3u dx 



6 * 2 (c6* + <$* 2 tana) 



(eqn 2.23) 



where a = fin half angle (radians) 
e = width of the trough (ft) 

Note that the quantity in parentheses is the cross sectional 
area of the film condensate in the trough (See Figure 4) . 
Taking the derivative of equation (2.23) with respect to x 
yields the rate of change of mass flo;\r in the trough for a 
given axial increment. 



dm total 

dx 



P f u) 2 r 

3y 



— [4^ (c<S* 3 + <5* 4 tana) ] 

dx dx 



(eqn 2.24) 



Equation (2.24) represents the rate of change of the 
total mass flow rate with respect to x in the x-direction. 
This equation must be coupled with the energy equations for 
the fin and trough to develop a representation of the film 
profile in the trough. 
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3. Mass Flow in the Z-Direction 



Examining an infinitesimal fluid element on the sur- 
face of a fin for any axial increment of width Ax, as shown 
in Figure 6, the momentum equation in the z-direction becomes: 




3P 
3 z 



p £ u) 2 rcosa 



(eqn 2.25) 



where r_ = shear stress in the z-direction (lbf/ft 2 ) 

Zf 

z = co-ordinate measuring distance along the surface 
of the fin (ft) 

Neglecting dP/dz based on assumption (b) , and integrating 
equation (2.25) from i, to 0 and y to 6 yields: 



x z = u 4- = p £ o) 2 rcosa(5-y) 



(eqn 2.26) 



where w = fluid velocity in the z-direction (ft/hr) 

6 = fin film thickness along the surface of the fin (ft) 
Note, 6 , the fin film thickness should not be confused with 
<5*, the film thickness in the trough. Integrating equation 
(2.26) from 0 to w and 0 to y provides the following expression 
for fluid velocity: 

p £ 6 2 rcosa 2 

w = (S(z)y - \ ) (eqn 2.27) 
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Figure 6. Cross Section of Infinitesimal Fluid Element on 
Fin Surface of Axially Finned Condenser 
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This relationship may be used to find the average fin fluid 
velocity w: 



• w = 




p^.a) 2 rco sa 5 ( z) 2 
3y 



(eqn 2.28) 



The mass flow rate in the z-direction along the surface of the 
fin for a given axial increment is given by: 

m^ n = P ^wdA (eqn (2.29) 

where dA = the cross sectional area of the fluid flowing along 
the fin surface (ft 2 ) 

Substituting equation (2.28) into equation (2.29) yields: 

P£ 2 w 2 rcosa S 3 (z)dx 

m fin = (eqn 2.30) 

3p 

This equation is identical to equation (IS) of Schafer's [Ref. 

3] analysis if the condenser cone half angle (0) is set equal 
to 0 which is the case for a cylindrical condenser. 

4 . Energy Equation for the Trough Condensate 

An energy balance on an infinitesimal fluid element in 
the trough of an axial increment of width dx with surface area 
e • dx yields the following expression for heat transfer by 
conduction: 

trough * kf(Tsat 1 Tw)£dX (e <t n 2 ' 31) 
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Note also, that the trough heat transfer rate is given by: 



^trough h fg ^trough 



(eqn 2.32) 



Combining equations (2.31) and (2.32) and dividing by dx re- 
sults in the following: 



^trough 

dx 



k f< T sat - 



fg 



(eqn 2.33) 



Equation (2.33) is an expression for incremental change in mass 
flow rate with respect to x due to condensation in the trough 
region . 

5 . Energy Equation for the Fin Condensate 

An energy balance on a differential element of surface 
area dx*dz yields the following relationship for differential 
heat into the fin: 



do r • h dm c • 
H fin fg fin 



k r (T *. 
f ^ sat 



- T fin (z))dxdz 

6(z) 



(eqn 2.34) 



where T^ n (z) = fin surface temperature at some position z 

along the surface of the fin (degrees F) 

Since the fin condensate mass is assumed to flow only in the 
z-direction, equation (2.30) is differentiated with respect 
to z and substituted into equation (2.34). After substitution 
and rearrangement, the following equation results: 
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5(z) 3 d<5 (z) 



(eqn 2.35) 



V T sat - T 

p f 2 “ 2r h fg 






cosa 



Applying assumption (c) i.e., T^ n (z) equals T avcr for 
all z and integrating equation (2.35) from 0 to 6 and 0 to z 
yields the following relationship for fin film thickness <5 (z) : 



6(z) = 



4 k r (T . - T ) y z l 

f ^ sat avg^ 

p^ 2 o) 2 r h^g cosa 



1 / "4 



(eqn 2.36) 



where T = average fin surface temperature (degrees F) . 

v o 

Substituting equation (2.36) into equation (2.30) and solving 
for rate of change of mass flow rate of the fin with respect 
to x for an increment of width dx yields : 



dm 



fin 



dx 



2p£ 2 w 2 rcosa 

3y 



4 k £( T sat 



T )pz* ■ 
avg - ^ M 



P£ 2 w 2 r h^ a cosa 



(eqn 2.37) 



where z* = z-6*/cos(a). Note, z* is the distance along the 
surface of the fin from the fin tip to the trough film thick- 
ness (6*). Note also, that the right hand side of equation 
(2.37) is multiplied by two; this accounts for mass flow from 
the fins on both sides of the trough. 

6 . Continuity Equation 

For any axial increment of length dx, continuity dic- 
tates that: 



39 



dm 



total 



dm 



fin 



dm 



dx 



dx 



trough 

dx 



( e qn 2.38) 



Substituting equations (2.24), (2.33) and (2.37) into equation 
(2.38) and rearranging yields: 



«* [jf* (£«* = * «*‘taiK»)l 



5 k d T sat ' 
p£ 2 w 2 r F fg 



- 26*cos 



•1 



4 k r (T . - T )pz* ■> 
f v sat avg ; 

p^ 2 w 2 r h£g cosa 



3/4 



(eqn 2.39) 



Equation (2.39) can be solved using the Finite Element Method 
formulation provided in Appendix A. The solution of this equa- 
tion provides the film thickness profile along the axial length 
of a cylindrical finned condenser. 

7 . Determination of the Heat Transfer Rate 

Once the film profile has been determined within the 
trough, the local convective heat transfer coefficient can be 
found for the trough using the following relationship: 



h(x) t 



rough 



5*(x) 



(eqn 2.40) 



In a similar manner, the local heat transfer coefficient along 
the surface of the fin can be found by: 



h(z) 



fin 6(z) 



(eqn 2.41) 
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The differential heat transfer rate for any fin section, as 
shown in Figure 4, of axial incremental length dx is: 
a) for the trough: 



. . CT sat - Tjsdx 

q trough h(x) trough 



(eqn 2.42) 



where £ # dx is the surface area of the trough, and 
b) for the fin surface: 



dq 



fin 





- T 
U) 



avg 

fin 



) dxdz 



(eqn 2.43) 



where z is the surface length of the fin. 

The total differential heat transfer rate per axial increment 
is found by summing equation (2.42) and (2.43) for the total 
number of fins. That is: 

NFIN 

dq total = 2 (dq fin * ^trough 0 (eqn 2 ' 44) 



where NFIN is the total number of axial fins . 

In a similar manner, the total heat transfer rate for the entire 
finned condenser can be found by the following relationship: 



NDI V 

^total ^ ^total 



(eqn 2.45) 



where NDIV is the total number of axial increments. 
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III. COMPUTER CODE DESCRIPTION 



A. GENERAL DESCRIPTION OF CODE 

The computer code consists of a main body and eight sub- 
routines. Basically, the code which is provided in Appendix 
C is a modification of Purnomo's [Ref. 5] code. The function 
of each subroutine used in the code is as follows: 

a) "CORRES" established the correspondence between the 
local and global nodal points used in the Finite Element 
Method solution for the two-dimensional steady state 
heat conduction problem. In so doing, "CORRES" also 
numbers all elements and nodal points in the finite ele- 
ment model and assigns local nodal points to each of the 
elements. In addition "CORRES" also defines major ele- 
ment numbers used in other subroutines as control 
parameters . 

b) "COORD" defines the x and y coordinates for all nodal 
points in the finite element heat conduction problem 
model. 

c) "DLSTAR" determines the film thickness (6*) on the sur- 
face of a smooth condenser or in the trough in a finned 
condenser . 

d) "HTCOEF" determines the heat transfer coefficient for 
all convective surface elements. 
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e) "FORMAF" formulates the Finite Element Method equations 
for the two-dimensional steady state heat conduction 
problem. 

f) "BANDEC" is an equation solver for a symmetric matrix 
which has been transformed into banded form. "BANDEC” 
will return the solution to the two-dimensional heat 
conduction problem. 

g) • "HTCALC" determines the elemental, incremental and total 

heat transfer rates. 

h) "DELCRV" determines the condensate film profile in a cyl- 
indrical condenser. 

Two additional Naval Postgraduate School computer library 
routines are also used in the code: 

a) "DPOLRT" is a nonIMSL double precision library routine 
that determines the roots of a real polynomial. This 
routine is called by "DLSTAR" to determine the film 
thickness for the succeeding increment in the analysis 
of a tapered condenser. 

b) "LEQT2F" is an IMSL double precision libary routine that 
solves a set of simultaneous linear equations. This 
routine is called by "DELCRV" to solve the Finite Element 
Method equations for the cylindrical condenser film pro- 
file problem. The resulting film profile is then used 

in the heat conduction analysis. 

In order to use the computer code to analyze heat transfer 
in a rotating heat pipe, nine data cards are required. A user’s 
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guide describing these data cards and required input is pro- 
vided in Appendix B. The input data, describing the geometric 
configuration of the rotating heat pipe as well as the operating 
parameters determines which solution technique is utilized in the 
analysis. The solution technique for each of the four con- 
denser geometries, i.e., tapered- smooth , tapered- axially finned, 
cylindrical - smooth and cylindrical - axially finned is different. 

In all cases however, the Finite Element Method is used to 
solve the two-dimensional steady state heat conduction problem. 
This solution is the one developed by Purnomo [Ref. 5] and 
has not been modified. Details of the development of this 
solution are described in detail in Purnomo ' s thesis [Ref. 5] 
and will not be repeated here. This being the case, each of 
the four solution techniques will now be discussed in detail. 

B. INTERNALLY FINNED TAPERED CONDENSER SOLUTION 

The complete development of this solution technique is 
provided in Purnomo ' s [Ref. 5] thesis and will not be rede- 
veloped. When an equation is required for clarity, the equa- 
tion in final form will be provided. Where there is a 
modification to Purnomo ' s [Ref. 5] code, this modification 
will be noted. 

The condenser of the internally finned- tapered condenser 
is divided into NDIV axial increments. These axial incre- 
ments are then subdivided circumferentially into ZFIN number 
of subincrements where ZFIN is the total number of fins. 



44 



These subincrements are then divided in half to form the 
basic symmetric unit, one-half of a fin-trough section as 
shown in Figure 1. This unit section is then divided into 
a number of linear triangular elements. The number of ele- 
ments depend on the input parameters. The only limitation 
is that the number of system nodal points must hot exceed 
100; otherwise, certain variables, e.g., x and y, would 
exceed allotted storage values. Figure 7 shows a unit 
section subdivided into 25 elements. After this unit is 
subdivided, each nodal point is assigned an x and y coordi- 
nated based on the geometric input parameters. 

To start the iteration, two initial values are required: 
1) an initial temperature for the nodal points along the 
internal convective boundary, and 2) an initial trough film 
thickness (6*). The initial temperature is provided as an 
input parameter and the initial trough film thickness at the 
first increment is provided by a relationship taken from an 
analysis by Sparrow and Gregg for condensation on a rotating 
disk [Ref. 12 ] . 

Once these values are known, the heat transfer coeffi- 
cient for the internal convective surface elements are found 
using the following relationships: 



1 / 



h(z) 



fin 



_ k f 



6Tz) 



kf Pf w 2 (r+xsin0) h£gCos0cosa 
4y (-AA- z J /3-BB- z 2 /2+(T sat -T 1 ) z 



(eqn 3.1) 
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a = 0.02953 inches 
b = 0.05906 inches 
c = 0.04412 inches 
t = 0.05118 inches 



Figure 7. Axially Finned Condenser Symmetric Section Sub- 
divided into 25 Linear Triangular Finite Elements 
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where x = distance from condenser end (x=0.0) to midpoint 
of increment (ft) 

= fin apex temperature (degrees F) 

AA, BB = constants in the parabolic temperature determined 
by a Langrangian fit. 

For the trough surface elements : 



t, r ^f 

trough ” 6 * (x) 



(eqn 3.2) 



These heat transfer coefficients, along with the thermal 
conductivity and x and y coordinates of the nodal points are 
used to form the Finite Element Method equations for the two- 
dimensional heat conduction problem. The equations are then 
solved to yield a temperature distribution in the symmetric 
section . 

The above iteration is repeated, where, now the solution 
temperature distribution from the previous iteration is used 
to calculate the heat transfer coefficients along the con- 
vective surface fin elements as well as a new 6* using Sparrow's 
and Gregg's relationship [Ref. 12], This neiv 6*, in turn, is 
used to determine the heat transfer coefficients of the trough 
elements . 

Again, the Finite Element Method subroutines will yield 
a temperature distribution for the symmetric section. At 
this point, the nodal point temperatures are checked for con- 
vergence using the following relationship: 
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< CRIT i = 1,2, . . . NSNP 



(eqn 3.3) 



T. • - T. • , 
i.J 1,3-1 

T- • 



where NSNP = number of system nodal points 
j = present iteration 
j -1 = previous iteration 
CRIT = convergence criterion 

If this convergence test is successful for all nodal points 
the increment is considered solved. If any one nodal point 
fails, the iterative process is repeated until convergence is 
met. The convergence test in equation (3.3) is different than 
the one used by Purnomo [Ref. 5] in his thesis. Purnomo com- 
pared incremental heat transfer rates per unit of condenser 
length, Q ^ , rather than temperature as is done in the modified 
code . 

If convergence is met, the heat transfer rate is determined 
From this heat transfer rate, the incremental mass flow rate is 
determined by the following equation: 



™total 



2 Qi Ax 



(eqn 3.4) 



where = heat transfer rate per unit length (Btu/hr-ft) 

Ax = incremental width (ft) 

Using this value of incremental mass flow rate determined 
by equation (3.4), the following equation is used to calculate 
the subsequent interval's trough condensate film thickness ( 5 *) 
with a polynomial rootfinder subroutine: 
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m 



P f 2 co 2 (r+xsin0) 6 * 2 (x)sin0(<5 * (x) e + 6* 2 (x) tana) 



total 



3]i 



(eqn 3 . 5) 



where e = trough width (ft) 

This resulting value of S*(x) is then defined as the trough 
film thickness for the next increment. In addition, the solu- 
tion temperature distribution from the previous iteration is 
used as the starting temperature distribution for the next 
increment . 

This iterative process at each increment is repeated until 
convergence is met, and, is continued at each increment until 
the entire length of the condenser has been transversed. In- 
cremental heat rates are then summed to yield the total heat 
transfer rate. That is: 

NDIV 

Qtotal = 2 * ZFIN * E Qi' Ax (eqn 3.6) 

i- 1 

where ZFIN is the total number of axial fins. 

Once the total heat transfer rate has been determined, the 
problem is solved and pertinent data is provided as output. 

C. SMOOTH TAPERED CONDENSER SOLUTION 

The heat pipe condenser is divided into NDIV number of 
axial increments as in the finned- tapered condenser solution. 
These axial increments are then subdivided into 360 segments 
of equal length; these segments are the basic symmetric unit. 
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This unit section, is divided into a number of linear tri- 
angular elements with the same limitation as before; the 
number of system nodal points must not exceed 100. The 
system nodal points are then assigned x and y coordinates 
based on the input geometric parameters . 

To start the iterative process, as in the finned- tapered 
case, the initial value of temperature which is an input 
parameter is used to solve for the 'initial value of fin 
thickness (6*) based on the Sparrow and Gregg analysis [Ref. 



Once this initial value of 5* is known, the heat transfer 
coefficients for the internal convective elements can be 
determined using equation (3.2). These heat transfer coef- 
ficients are used in the Finite Element Method equations. 

The equations are solved yielding a temperature distribution. 
The iteration is repeated until convergence is met, just as 
in the finned- tapered case. 

When convergence is met, that is equation (3.3) has been 
satisfied, a new film thickness 5*(x) can be found by one of 
the following equations for 6*(x): 



12 ]. 




1=0 (eqn 3.7) 



or 




1 = 0 



(eqn 3.8) 
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<5*U) = 



3k gP (Tgat ~ T w ^ 
2p f 2 w 2 rsin 2 0h f 



[1 



{- 



(r+xs in0) 
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iA 



(eqn 3.9) 



where Sh = 



P £ 2 (w 2 r-g)sin0Pi £g x 



4 U ^ff T S at - V 



Dr = 



Pj-T hj- X 2 COS0 
f V fg 

u k f( T sat - T w> 



P^-V COS0 
Rev = — 

y 

g = acceleration due to gravity (ft/hr 2 ) 
v = vapor velocity (ft/hr) 

T = local wall temperature (degrees FO 
w 

t = shear stress vapor-liquid interface (lbf/ft 2 ) 
Equation (3.7) defines the film thickness distribution for 
a smooth tapered rotating heat pipe derived by Daniels and 
Al-Jumaily [Ref. 13]. This equation takes into account the 
drag effects of counter- flowing vapor. Equation (3.8) is a 
modification of Equation (3.7) neglecting the drag losses. 
Equation (3.9) was developed by Ballback [Ref. 1]. This equa- 
tion also neglects drag. 

Depending on a particular control parameter which is part 
of the input data (See Appendix B) one of these equations is 
used to solve for the film thickness (6*(x)) for the next 
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increment. In addition, the solution temperature distribution 
from the previous iteration is used as the starting temperature 
distribution for the next increment. 

This iterative process at each increment is repeated until 
convergence is met, and is continued at each increment until 
the entire length of the condenser has been transversed, just 
as in the finned- tapered case. Total heat transfer rates are 
then determined by summing the incremental heat transfer rates 
for the entire length of the condenser by the following 
relationship : 

NDIV 

Q t otal = 36 °* Qi Ax Oqn 3.10) 

D. SMOOTH CYLINDRICAL CONDENSER SOLUTION 

As in the smooth- tapered case, the condenser is first 
divided axially, then it is divided circumferentially into 
360 segments of equal length. These segments are the basic 
symmetric unit section to be considered. Again the symmetric 
unit is subdivided into linear triangular elements and x and 
y coordinates are assigned to the system nodal points. 

To begin the iteration, an initial temperature estimate 
which is an input parameter is assigned to the convective sur- 
face nodal points. Using this initial temperature estimate, 
the maximum film thickness, 5* „ which is located at x = 0 , 
is determined. This maximum film thickness value is then 
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used as one of the boundary conditions in the solution of 
equation (2.18) using the Finite Element Method. Equation 
(2.18) is repeated here for reference. 




* d T d * * 3 



_ _ 5kf (T sa t ~ h 

p f 2 w 2 r h f0 

-L & 



(eqn 3.11) 



The finite element solution of Equation (3.11) provides the 
film thickness profile (S*(x)) at the midpoint of each incre- 
ment along the length of the condenser, but only the first 
increment value is applied at this point of the analysis. 

Once 6*(x) is known at the first increment the heat transfer 
coefficients for the internal convective surface elements can 
be determined using equation (3.2). The steady state heat con- 
duction problem is then solved and a temperature distribution 
for the unit section results. This process is now repeated. 

A new maximum film thickness, film profile and thus <5*(x) at 
the first increment is found based on the solution temperature 
distribution from the first iteration. With this new value for 
6*(x), the heat conduction problem is again solved for a new 
temperature distribution. This iterative process, involving 
solution of the film profile with each iteration is continued 
until temperature convergence is met at the first increment. 
When convergence is met, the film profile determined on the 
iteration in which convergence was met is used to provide the 
values of 6*(x) for the remaining increments along the length 
of the condenser. 
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Using this predetermined value of 5*(x) , the iterative pro- 
cess is continued at each increment until temperature conver- 
gence is reached. As convergence is reached at each increment, 
the internal wall temperature, which is the same for all inter- 
nal wall nodal points of the section, is stored for future 
application. 

When convergence is reached at the final increment, equation 
(3.11) is once again solved for the film profile. It should be 
noted, however, that the right hand side of equation (3.11) is 
temperature dependent. It should also be noted that the temper- 
ature varies axially along the length of the condenser. In 
order to account for this temperature dependence and temperature 
variation, the right side of equation (3.11) is now determined 
for each increment using the wall temperatures that were stored 
at each increment. The finite element solution of the film pro- 
file will now account for the temperature variation along the 
length of the condenser. 

This final film profile provides the value of 6*(x) for each 
increment. The iterative process of solving for heat transfer 
coefficient, temperature distribution, and temperature conver- 
gence is continued for each increment until the final increment 
is reached. 

Equation (2.13) provides a relationship for total mass flow 
rate as a function of position. The total mass flow rate at 
the overfall into the evaporator is given by: 
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(eqn 3.12) 



2Ttpf 2 co 2 r 



m 



d<S*(L) 

dx 



* 3 

3 CL) 



where 5*(L) and d<$*(L)/dx are the values for the film thickness 
and rate of change of film thickness with respect to x at the 
evaporator end of the condenser. Another relationship for mass 
flow rate based on the steady state heat conduction solution 
is given by: 



m = 



NDIV 0-Ax 
360 * 2 

i " 1 h fg 



(eqn 3.13) 



If, in fact, the solution of equations (3.12) and (3.13) are 
equal, then the mass flow rate of the condensate returning to 
the evaporator is equal to the mass flow rate of the vapor 
being condensed on the surface of the condenser. Or, to put 
it another way, continuity is satisfied. 

It should be noted that the film profile maintains the same 
basic shape, that is, o * ( x ) at x=0 is always equal to 5* ax and 
decreases to a specific minimum value at x=L. This being the 
case, if the maximum film thickness is varied, the film thick- 
nexx profile will vary in the same manner. For example, if 
the maximum film thickness is increased, the entire profile 
will also increase. This will result in an increased internal 
thermal resistance and thus lower heat transfer rate. As a 
result of the lower heat transfer rate, the mass flow rate as 
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determined by equation (3.13) will be less. At the same time, 
however, the greater film profile will result in a greater 
value of S*(x) at the overfall, x = L. Yet, since the profile 
maintained the same basic shape, the derivative at the overfall 
remains relatively constant. Thus the mass flow rate as deter- 
mined by equation (3.12) will increase. A decrease in the 
maximum film thickness will result in an opposite effect to 
the mass flow rates determined by equations (3.12) and (3.13). 

This being the case, the mass flow rates, as determined 
by equations (3.12) and (3.13) are now compared to determine 
if the film profile is in fact the solution profile to the 
problem. If continuity is not satisfied, the maximum film thick- 
ness is varied and the entire iterative process is restarted. 

This process is continued until the film profile mass flow rate, 
equation (3.12) converges towards the heat transfer mass flow 
rate, equation (3.13). When the absolute difference between 
these mass flow rates is less than a specific value, the re- 
sulting heat transfer rate is considered the solution to the 
problem. 

E. FINNED CYLINDRICAL CONDENSED SOLUTION 

As in the finned- tapered case, the condenser is divided 
axially and circumferentially into the basic symmetric section 
as shown in Figure 1. This unit is then subdivided into linear 
triangular elements and x and y coordinates are assigned to 
each nodal point. 
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To begin the iterative process, as before, an initial temp- 
erature estimate is assigned to each nodal point along the in- 



ternal convective surface. An initial trough film profile is 
determined, using equation (3.11) and this initial temperature 
estimate. In this case, however, the maximum film thickness 
(6* ax ) is not calculated but is an input parameter. 

Once o*(x) is known at the first increment, the internal 
heat transfer coefficients are determined, using equations (3.1) 
and (3.2). These values are then used in the Finite Element 
Method solution of the steady state heat conduction problem. 

A temperature distribution is determined and the iteration is 
repeated until temperature convergence is met. Note that a 
new film profile is determined for each iteration. 

As in the cylindrical-smooth condenser case, once conver- 
gence is met at the first increment, the film profile that was 
determined for the iteration prior to convergence at the first 
increment is then used to provide the film thickness 6*(x) for 
the remaining increments. 

At the final increment, a new film profile for a finned 
cylindrical condenser is then determined by solving the follow- 
ing equation developed in Chapter II: 



* d r d-5 
3 dx L dx 



U* 3 e + 5 * 4 tana) ] = - 5k f( T sat ~ Tw)ye 

P f 2 u 2 r h fg 




(eqn 3.14) 
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A description of the solution of this differential equation 
using the Finite Element Method is provided in Appenxix A. 

The iterative process of finding the solution temperature 
distribution for all increments is then repeated until the 
length of the condenser has been transversed. 

At this point, the total mass flow rate of the condenser 
returning to the evaporator is determined by the following 
relationship : 

p f 2 w 2 r 

m = jj-j — ■■ (e 6 * 3 + 6 * 4 tana) * ZFIN (eqn 3.15) 



where ZFIN is the number of axial fins. 

I 

This mass flow rate is compared to the mass flow rate given 
by the following equation: 



m 



ZFIN*2 



NDI V 
* Z 
i = l 




(eqn 3.16) 



Just as in the smooth-cylindrical condenser case, if the 
absolute difference between the two mass flow rates is less 
than a mass flow convergence criterion, the problem is con- 
sidered solved. If not, „ is varied and the entire itera- 
tive process is started again. As in the smooth- cylindrical 
condenser, varying 5 * will have the same effect on the film 
profile and the heat transfer rate. Since the temperature dis- 
tribution along the condenser has been solved once and closely 
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approximates the final solution to the problem, on successive 
iterations, equation (3.14) may be used to solve for the film 
profile rather than equation (3.11). Equation (3.11) was used 
on the first iteration because the finite element solution con- 
verges more quickly than equation (3.14) when an estimated 
temperature is used. 

A word of caution is required. The solution of equation 
(3.14) is highly sensitive to the value of 6* ax « If the initial 
value of 6* is inconsistent with the actual solution, e.g. 
too small, the Finite Element Method solution of the film pro- 
file will not converge. If this is in fact the case, the prob- 
lem will be automatically terminated. A new value of 6* should 
then be chosen and the problem restarted. 

One additional topic which should be addressed is the rec- 
tangular fin solution process. The rectangular fin profile is 
a slight modification of the finned- tapered or finned-cylindrical 
condenser solution. The only variation is that the top elements 
of the rectangular fin are assigned heat transfer coefficients 
of 0.0. Thus, the tip of the rectangular fin is considered 
adiabatic and no heat is transferred through this top face. 

Other than this modification, the solution techniques are the 
same as for the finned cases addressed above. 
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IV. RESULTS AND DISCUSSION 



Prior to the heat transfer analysis of the various con- 
denser configurations, it was necessary to verify the finite 
element solution of the film profile. The development of this 
solution is discussed in Appendix A. 

Leppert and Nimmo [Refs. 8 and 9] had developed an analy- 
tical solution for film condensation on a horizontal plate at 
a constant surface temperature. Their analysis and resulting 
differential equation is identical to the development in Chapter 
II for a smooth cylindrical condenser if the acceleration due 
to gravity is replaced by a radial acceleration term. This 
being the case, this modification was made and the analytical 
solution was used as a reference for comparison. 

For the test runs of the finite element solution, a constant 

surface temperature was assumed. In addition, a value for the 

maximum film thickness (<$* ax ) and minimum film thickness (o* an ) 

at the overfall were required. The value for < 5 * was determined 

max 

based on a relationship developed by Leppert and Nimmo [Ref. 8 ]. 
The value for S m ^ n was arbitrarily chosen. 

For identical geometry, surface temperature and maximum and 
minimum film thicknesses, the results of both analyses were 
identical. In order to develop confidence in the finite element 
solution, S* in was varied from 0.10 • 6 * ax to 0.97 • 6 * ax » The 
resulting profiles agreed at all locations along the length of 
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the condenser. However, when the temperature was permitted to 
vary along the length of the condenser in the finite element 
solution and a corresponding average temperature was used in the 
analytical solution of Leppert and Nimmo [Refs. 8 and 9], the 
profiles were no longer in agreement. This was to be expected, 
particularly in the case where there was a sixteen degree 
Fahrenheit variation along the length of the condenser. This 
substantial temperature variation resulted in a significant 
variation of fluid properties which would account for a differ- 
ence in film profile. 

Due to the agreement between the finite element solution 
and the analytical solution of Leppert and Nimmo [Refs. 8 and 

V 

9] for a constant surface temperature, it was decided that the 
finite element solution does provide a satisfactory representa- 
tion of the film profile. This being the case, the finite 
element solution was then incorporated into the code to provide 
the film profile for the cylindrical condenser. 

Once the finite element solution of the film profile was 
verified, the heat transfer analysis could be accomplished. 

The analysis considered both copper and stainless steel con- 
densers with the following four geometries: a) tapered- smooth , 

b) tapered-axial ly finned, c) cylindrical-smooth and d) cyl- 
indrical-axially finned. Table I lists the geometric para- 
meters held constant for all analyses. In all cases, the 
working fluid was water. 
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TABLE I 



Condenser Geometric Parameters 
Held Constant During All Analyses 

condenser length = 8.500 inches 

minimum radius = 0.51575 inches 

wall thickness = 0.05118 inches 

In addition, the following geometric parameters were also 
utilized when required. This requirement was based on the 
condenser geometry being considered, i.e., tapered-axially 
finned. 



TABLE II 

Condenser Geometric Parameters 
Applied as Required 



height of fin 


= 


0.05906 


inches 


fin half angle 


= 


26 . 565 


degrees 


condenser cone 








half angle 


= 


1.00 


degrees 



In the analysis, the heat transfer rate was determined for 
the four different geometries listed above for both copper and 
stainless steel. The ambient temperature was set at 60.0°F 
and the heat transfer rate was determined for each possible 
combination of the operating parameters given in Table III. 
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TABLE III 



Operating Parameter Matrix 



Rotat ional 



Speed 

(RPM) 



(btu/hr- ft 2 -F) 



Heat Transfer 
Coefficient 



Saturation 
Temperature 
(degree F) 



700.0 



100 . 0 



90.0 



1400.0 



500.0 



120.0 



2800.0 



1000.0 



150.0 



180 . 0 



Thus, for each condenser geometry, there was a total of 72 
analyses, 36 for the case of the condenser with a copper wall 
and 36 for the case of the condenser with a stainless steel 
wall . 

The first condenser geometry considered was a smooth con- 
denser. Figures 8-13 compare the heat transfer rates of smooth 
cylindrical condensers with those of smooth tapered condenser. 

In particular. Figures 8,9, and 10 indicate the results of the 
analyses of smooth copper condensers at rotational speeds of 
700, 1400 and 2800 revolutions per minute (RPM) respectively. 
Figures 11, 12, and 13 are for smooth stainless steel conden- 
sers at 700, 1400, and 2800 RPM respectively. For both stain- 
less steel and copper smooth condensers, the following general 
observations apply: a) For the same external heat transfer 

coefficient, the heat transfer rate for the cylindrical smooth 
condenser is less than the equivalent tapered condenser, b) As 
the external heat transfer coefficient increases, this difference 
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in heat transfer rate also increases. c) The rotational speed 
has a greater effect on the tapered condenser heat transfer 
rate. For example, the maximum heat transfer rate of a copper 
tapered condenser will increase by a factor of 1.67 when the 
rotational speed is increased from 700 to 2800 RPM. In the 
cylindrical copper condenser, for the same change in rotational 
speed, the heat transfer rate only increases by a factor of 1.51. 

These same observations hold true for the smooth stainless 
steel condensers. But, due to a greater thermal resistance in 
the wall of the stainless steel condenser, the heat transfer 
rates are less for all cases considered. It should be noted 
that the thermal conductivity of stainless steel is only 4% the 
thermal conductivity of copper. This accounts for the increased 
thermal resistance. 

Figures 14, 15, and 16 compare axially finned cylindrical 
with axially finned tapered copper condensers at 700, 1400 and 
2800 RPM respectively. Note that the heat transfer rates of 
the cylindrical condensers are only slightly less than those of 
the tapered condensers. This is because the heat is primarily 
transferred through the extended surface, i.e., the fin. Thus 
the film condensate in the trough has less effect on the heat 
transfer rate. In fact, the average difference in heat transfer 
rate for 700 RPM and an external heat transfer coefficient of 
100 Btu/hr-ft 2 -F is 12.65%. This difference increases to 15.6% 
as the heat transfer coefficient is increased to 1000 Btu/hr- 
ft 2 -F. However, as the rotational speed is increased to 2800 
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RPM, the corresponding average differences in heat transfer 
rates decrease to 12.58% and 13.33% respectively. 

Figures 17, 18, and 19 compare the heat transfer rates of 
axially finned stailness steel cylindrical condensers with 
those of axially finned stainless steel tapered condensers 
at 700, 1400 and 2800 RPM respectively. Note the difference 
in heat transfer rate increases as the external heat transfer 
coefficient increases more than in the case of the copper 
condensers. At the low heat transfer coefficient, the limit- 
ing thermal resistance is that of the external surface (1/h . 

However, as the external heat transfer coefficient is increased, 
the limiting thermal resistance becomes that of the wall due 
to the low thermal conductivity of stainless steel. 

Another observation to be noted is the fact that the rota- 
tional speed has a greater effect on the heat transfer rate of 
the cylindrical condensers than on the tapered condensers . As 
the rotational speed increases, the film thickness decreases 
due to the greater centrifugal force exerted on the film. For 
the cylindrical condenser, the maximum film thickness is much 
greater than for a tapered condenser. For example, an axially 
finned stainless steel cylindrical condenser rotating at 1400 
RPM with an external heat transfer coefficient of 1000 Btu/hr- 
ft 2 -F'has a maximum film thickness twice that of a corresponding 
tapered condenser. This being the case, higher rotational 
speeds will have a greater effect on the greater film thickness 
and the difference in heat transfer rates will decrease. 
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Figures 20, 21, and 22 compare the heat transfer rates of 
copper smooth cylindrical condensers with copper axially finned 
cylindrical condensers at 700, 1400 and 2800 RPM respectively. 
The figures indicate that for low external heat transfer coef- 
ficients, little is to be gained by the addition of axial fins. 
But, as the external heat transfer coefficient increases, the 
advantage becomes significant. As an example, consider Figure 
21. For h=100 Btu/hr- ft 2 -F , axial finning increase the heat 
transfer rate by 16%. But, for an external heat transfer coef- 
ficient of 1000 Btu/hr-ft 2 -F, the heat transfer rate increases 
by 194%. 

Note also, that as rotational speed increases, the advan- 
tage to be gained by axial finning decreases. This can be 
explained by the fact that for an axially finned condenser, 
the majority of the heat is transferred by the fin surface. 

As the rotational speed increases, the film thickness in the 
trough decreases. This in turn will expose slightly more fin 
surface area. On the other hand, for the smooth condenser, 
the decrease in the film thickness will have a greater effect 
on heat transfer rate in that the thermal resistance of the 
film has decreased. Comparing the two geometries, the change 
in overall thermal resistances for the smooth condenser will 
be greater than that for the axially finned condenser account- 
ing for the slight decrease in advantage with increasing 
rotational speed. 
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Figures 23, 24, and 25 correspond to Figures 20, 21, and 
22 but for stainless steel condensers. The results are similar 
to the copper situation discussed above, by the effect caused 
by the increasing external heat transfer coefficient is not 
as dominant due to the high thermal resistance of the stainless 
steel wall material. 

Figure 26 compares a smooth cylindrical copper condenser 
with a smooth cylindrical stainless steel condenser at 1400 
RPM. As to be expected, the heat transfer rate of the copper 
condenser is greater than that of the stainless steel condenser 
due to the difference in the thermal conductivity of the two 
materials. The difference in the heat transfer rate is least 
for a low external heat transfer coefficient where the external 
thermal resistance is dominant. As the heat transfer coeffi- 
cient increases, the thermal resistance of the wall of the 
condenser becomes more important resulting in an increasing 
difference between the two condensers. 

Figure 27 provides a comparison of axially finned cylindrical 
copper and stainless steel condensers at 1400 RPM. The advantage 
of copper over stainless steel is obvious. Note that the heat 
transfer rate for the copper condenser at 500 Btu/hr-ft 2 -F is 
nearly identical to the stainless steel condenser heat transfer 
rate at 1000 Btu/hr/ - f t 2 - F indicating the advantage of copper 
over stainless steel. 

The final analysis compared the heat transfer rate of an 
axially finned copper condenser with a triangular fin profile 



67 



to the heat transfer rate of an axially finned condenser with 
a rectangular fin profile. The rectangular fin was assumed to 
have an adiabatic tip. This comparison was accomplished for 
both cylindrical and tapered condensers. Table IV lists the 
parameters used in the analysis. These parameters are in addi- 
tion to those listed in Table I. This comparison was only 
accomplished for the one set of operating parameters listed in 
Table IV. Note that the operating parameters chosen were the 
median values . 



TABLE IV 

List of Parameters Used in 
Rectangular/Triangular Fin Profile 



Heat Transfer Coefficient 
Rotational Speed 
Saturation Temperature 
Fin Height 

Table V lists the results 



Comparison 

500 Btu/hr-ft 2 -F 
1400 RPM 
120.0 degrees F 
0.05906 inches 
of the analyses . 



TABLE V 

Results of Rectangular/Triangular 
Fin Profile Comparison 



CONDENSER GEOMETRY 


FIN PROFILE 


Q (Btu/hr) 


tapered 


triangle 


6168.93 




rectangle 


6143.03 


cylindrical 


triangle 


5358.4 




rectangle 


5338.00 
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Note that the heat transfer rate varied by only 0.4% for both 
the tapered and cylindrical condensers. The reason for this 
insignificant variation is provided in Table VI and VII. 

Table VI lists the convective surface temperature distribution 
of a symmetric unit section of a tapered condenser at the middle 
increment of the condenser. Temperature location 1 is located 
at the tip of the fin. For the rectangular profile, temperature 
location 1 is located at the intersection of the adiabatic 
surface (the tip of the fin) , and the vertical surface of the 
fin. Temperature location 5 is located at the base of the fin. 
Locations 6-9 are located in the trough region and the remaining 
temperature locations are along the external surface of the 
section. Thus temperatures 1-5 provide the temperature distri- 
bution along the convective surface of the fin. Note that the 
temperature distribution is lower for the rectangular fin. 

Thus, the driving force for heat transfer, the temperature 
difference between the saturation temperature and the surface 
temperature is greater for the rectangular fin. Thus, in spite 
of the fact that the fin surface area for heat transfer has de- 
creased by 11%, the average surface temperature difference has 
increased by 25%. It should also be noted that the average 
heat transfer coefficient for the fin surface for the rectan- 
gular fin also decreases by 12%. Hoxvever, the increase in temper- 
ature difference is the dominant change that only allows a 0.4% 
decrease in heat transfer rate. 
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TABLE VI 



Surface Temperature Distribution for Rectangular 
and Triangular Axially Finned Tapered Copper 
Condensers at Middle Increment of Condensers 



Temperature Triangular Fin Rectangular Fin 

Location Temperature (F) Temperature (F) 



1 


118.85 


117.95 


2 


118.45 


117.88 


3 


118.04 


117.68 


4 


117.60 


117.36 


5 


117.05 


116.59 


6 


116 . 82 


116.45 


7 


116.61 


116.35 


8 


116.59 


116.35 


9 


116.48 


116.54 


10 


116.47 


116.53 


11 


116 . 46 


116.23 


12 


116.46 


116.21 


13 


116.44 


116.20 


14 


116.42 


116.17 


15 


116.31 


116.12 


16 


116.31 


116.07 


17 


116.27 


116.04 


18 


116.26 


116.03 
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TABLE VII 



Surface Temperature Distribution for Rectangular 
and Triangular Axially Finned Cylindrical Copper 
Condensers at Middle Increment of Condensers 



Temperature 

Location 


Triangular Fin 
Temperature (F) 


Rectangular 

Temperature 


1 


118.55 


117.55 


2 


118 . 09 


117.44 


3 


117.62 


117.25 


4 


117.10 


116.87 


5 


116.55 


116 . 36 


6 


116 . 38 


116.18 


7 


116.28 


116.07 


8 


116.22 


116.01 


9 


116.20 


115.99 


10 


116.03 


115.80 


11 


116 . 02 


115 . 79 


12 


116.01 


115 . 79 


13 


115.99 


115.77 


14 


115.97 


115 .75 


15 


115.94 


115.72 


16 


115.91 


115.67 


17 


115.88 


115.67 


18 


115.88 


115.66 
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Table VII indicates a similar situation exists in the 
cylindrical condenser. Again note the decrease in the average 
surface temperature of the fin for the rectangular profile. 
The surface area and heat transfer coefficient has decreased, 
but the increased temperature difference compensates for these 
changes, limiting the overall decrease in heat transfer rate. 
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Heat Transfer Rate (Btu/hr) 




Saturation Temperature (°F) 



1. h ext = 100 Btu/hr-ft 2 -F 

2. h ext = 500 Btu/hr-ft 2 -F 

3. h « 1000 Btu/hr-ft 2 -F 

C A L 



Figure 8. Heat Transfer Rate versus Saturation Temperature 
for Smooth Cylindrical and Tapered Copper Con- 
densers at 700 RPM 
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Heat Transfer Rate (Btu/hr) 




Saturation Temperature (°F) 



1. h 0xt = 100 Btu/hr- ft 2 -°F 

2. h gxt = 500 Btu/hr- ft 2 -°F 

3. h ext = 1000 Btu/hr -ft 2 -°F 

Figure 9. Heat Transfer Rate versus Saturation Temperature 
for Smooth Cylindrical and Tapered Copper Con- 
densers at 1400 RPM 
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Heat Transfer Rate (Btu/hr) 



16000 r 




Saturation Temperature (°F) 



1. h ext = 100 Btu/hr- ft 2 -°F 
2 - h 0Xt = 500 Btu/hr- ft 2 -°F 
3. h ext = 1000 Btu/hr- ft 2 -°F 

Figure 10. Heat Transfer Rate versus Saturation Temperature 
for Smooth Cylindrical and Tapered Copper Con- 
densers at 2800 RPM 
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Figure 11. Heat Transfer Rate versus Saturation Temperature 
for Smooth Cylindrical and Tapered Stainless 
Steel Condensers at 700 RPM 
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Heat Transfer Rate (Btu/hr) 




Saturation Temperature (°F) 



1. h ext = 100 3tu/hr-ft 2 -F 

2. ^ ext = 500 Btu/hr-ft 2 -F 

3. h ext = 1000 Btu/hr- ft 2 -F 



Figure 12. Heat Transfer Rate versus Saturation Temperature 
for Smooth Cylindrical and Tapered Stainless 
Steel Condensers at 1400 RPM 
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Heat Transfer Rate (Btu/hr) 




1. h ext = 100 Btu/hr-ft 2 -F 

2. h ext = 500 Btu/hr- ft 2 - F 

3. h t = 1000 Btu/hr-ft 2 -F 

GX L 



Figure 13. Heat Transfer Rate versus Saturation Temperature 
for Smooth Tapered and Cylindrical Stainless 
Steel Condensers at 2800 RPM 
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Heat Transfer Rate (litu/hr) 




Saturation Temperature (°F) 

1. h ext = 100 Btu/hr- ft 2 -°F 

2. ^ ext = 500 Btu/hr- ft 2 -°F 

3. h ext = 1000 3tu/hr-ft 2 -°F 



Figure 14. Heat Transfer Rate versus Saturation Temperature 
for Axially Finned Cylindrical and Tapered Copper 
Condensers at 700 RPM 
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Heat Transfer Rate (Btu/hr) 




1 

2 

3 



h ext = 100 Btu/hr-ft 2 -F 
h gx t = 500 Btu/hr-ft 2 -F 
h ext = 1000 Btu/hr-ft 2 -F 



Figure 15. 



Heat Transfer Rate versus Saturation Temperature 
for Axially Finned Cylindrical and Tapered Copper 
Condensers at 1400 RPM 
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Heat Transfer Rate (Btu/hr) 




Saturation Temperature (°F) 

1. h ext = 100 Btu/hr- ft 2 -°F 

2. h ext ~= 500 Btu/hr- ft 2 -°F 

3. h ext = 1000 Btu/hr- f t 2 - °F 

Figure 16. Heat Transfer Rate versus Saturation Temperature 
for Axially Finned Cylindrical and Tapered Copper 
Condensers at 2800 RPM 
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Figure 17. Heat Transfer Rate versus Saturation Temperature 
for Axially Finned Cylindrical and Tapered Stain- 
less Steel Condensers at 700 RPM 
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Heat Transfer Rate (Btu/hr) 




Saturation Temperature (°F) 

1. h gxt . = 100 Btu/hr- ft 2 -°F 

2. h ext = 500 Btu/hr- ft 2 -°F 

3. h gxt = 1000 Btu/hr- ft 2 -°F 

Figure 18. Heat Transfer Rate versus Saturation Temperature 
for Axially Finned Cylindrical and Tapered Stain- 
less Steel Condensers at 1400 RPM 
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Heat Transfer Rate (Btu/hr) 



Figure 19. 
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Heat Transfer Rate versus Saturation Temperature 
for Axially Finned Cylindrical and Tapered Stain- 
less Steel Condensers at 2800 RPM 
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Heat Transfer Rate (Btu/hr) 



Figure 2 
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0. Heat Transfer Rate versus Saturation Temperature 
for Smooth and Axially Finned Copper Cylindrical 
Condensers at 700 RPM 
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Heat Transfer Rate (Btu/hr) 




1. h ext = 100 Btu/hr- ft 2 -°F 

2. h ext = 500 Btu/hr-ft 2 -°F 

3. h ext = 1000 Btu/hr- ft 2 -°F 



Figure 21. Heat Transfer Rate versus Saturation Temperature 
for Smooth and Axially Finned Copper Cylindrical 
Condensers at 1400 RPM 
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Heat Transfer Rate (Btu/hr) 




Saturation Temperature (°F) 

1. h ext = 100 Btu/hr- ft 2 -°F 

2. h ext = 500 Btu/hr- ft 2 -°F 

3. h 0xt = 1000 Btu/hr- ft 2 -°F 

Figure 22. Heat Transfer Rate versus Saturation Temperature 
for Smooth and Axially Finned Copper Cylindrical 
Condensers at 2800 RPM 
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Heat Transfer Rate (Btu/hr) 




Saturation Temperature (°F) 



1. h ex t = 100 Btu/hr- ft 2 -°F 

2. h ext = 500 Btu/hr- ft 2 -°F 

3. h ext = 1000 Btu/hr- ft 2 -°F 



Figure 23. Heat Transfer Rate versus Saturation Temperature 
for Smooth and Axially Finned Stainless Steel 
Cylindrical Condensers at 700 RPM 
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Figure 24. Heat Transfer Rate versus Saturation Temperature 
for Smooth and Axially Finned Stainless Steel 
Cylindrical Condensers at 1400 RPM 
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Heat Transfer Rate (Btu/hr) 




1. h ext = 100 Btu/hr- ft 2 -°F 

2. h ext = 500 Btu/hr- ft 2 -°F 

3. h ext = 1000 Btu/hr- ft 2 -°F 



Figure 25. Heat Transfer Rate versus Saturation Temperature 
for Smooth and Axially Finned Stainless Steel 
Cylindrical Condensers at 2800 RPM 
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Figure 26. Heat Transfer Rate versus Saturation Temperature 
for Copper and Stainless Steel Smooth Cylindrical 
Condensers at 1400 RPM 
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Figure 27. Heat Transfer Rate versus Saturation Temperature 
for Copper and Stainless Steel Axially Finned 
Cylindrical Condensers at 1400 RPM 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The heat transfer rate of a cylindrical condenser is less 
than an equivalent tapered condenser. This decrease in heat 
transfer rate is most significant in a smooth condenser, where, 
depending on the external heat transfer coefficient and rota- 
tional speed, can be as great as 45%. This decrease in heat 
transfer rate becomes less significant for an axially finned 
condenser where the average decrease, fcr the range of heat 
transfer coefficients examined, was 13% for a copper condenser 
and 18% for a stainless steel condenser. When such factors 
as the cost and difficulty in manufacturing tapered axially 
finned condensers are considered, this 13% decrease in heat 
transfer rate becomes tolerable. From a practical standpoint, 
development and analysis of cylindrical axially finned con- 
densers should be encouraged by these results. If environ- 
mental conditions permit, copper should be preferred over 
stainless steel due to its exceptionally high thermal con- 
ductivity and resulting higher heat transfer rate. 

B. RECOMMENDATIONS 

1) Build and experimentally test both smooth and axially 
finned cylindrical condensers to obtain experimental data 
for comparison with results of this analysis. 
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2) Develop models for rectangular and trapezoidal fin 
profiles with non- adiabatic tips and incorporate into code 

3) Experimentally test axially finned cylindrical con 
densers with rectangular and trapezoidal fin profiles to 
obtain data for comparison with theoretical results. 
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APPENDIX A 



FILM PROFILE FINITE ELEMENT SOLUTION 
A. SMOOTH CONDENSER 



The analysis of the film profile in a smooth cylindrical 
condenser developed in Chapter II resulted in the following 
ordinary, nonlinear, second order, differential equation: 



6 * 



_d_ 

dx 



r d6* 

L clx“ 




3kf (Tsat 

p^ 2 'w 2 r K 



fg 



Tw) H 



(eqn A . 1) 



This equation can be rearranged and expanded to yield: 

6* 1 * - - + 36 * 3 (^— ) = -K (eqn A. 2) 

dx 2 dx 

where K = 3k f (T sat - T u )p 

p f ! u) 2 r lT fg 



The statement of the problem for the formulation of the 
Finite Element Method is: 



a** *±Al + 3 5 *> (S5H - -k 



dx' 



•d<S * 

•dx 



with the following boundary conditions 

a) at x = 0, = g * 

0 °max 

b) at x = 0 , d,s*/dx = 0 



(eqn A. 3) 
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The domain (length of the condenser) is divided into ele- 
ments of length Ax with the exception of the first and final 
elements which have length Ax/2, where Ax is the length of the 
condenser divided by the number of axial increments (NDIV) . 

A system nodal point is located at each end of an element. 
Thus, system nodal point 1 is located at x=0 and the last sys- 
tem nodal point, which is equal to the number of elements plus 
1 is located at the overfall into the evaporator, x=L. All 
internal nodal points are located at a position corresponding 
to the midpoint of the axial increment. 

Define the approximate value of the film thickness in the 
following manner: 



where "S’ = the approximate value of the film thickness (6*). 
= the global basis functions, 
n = the number of system nodal points, 
d = the solution vector. 

On an element level, equation (A. 4) becomes: 




h 




(eqn A. 4) 




4 

2 

i=l 




e 



(eqn A. 5) 



where g^ = the local basis functions. 
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Define two degrees of freedom at each nodal point, i.e., 
both ~6 and dS/dx are continuous. Thus, the local basis func- 
tions used in the finite element solution are: 



g, - 1 - ■* + ^ 

1 £ 2 £ 3 


(eqn A. 6) 


2C 2 £ 3 

g 2 - 5 - + , 

z l i z 


(eqn A. 7) 


_ _ 3? 2 2? 3 

§3 S, 2 S, 3 


(eqn A. 8) 


t 2 t 3 

= • r + f 


(eqn A. 9) 



where l = the length of an element. 

g^ and gj = magnitude basis functions. 
g 2 and g^ = slope basis functions. 

Substituting the approximate film thickness o" e into the 
differential equation (A. 3) results in: 



- 1 * d 2 o 0 . 3 ^d5 e 



dx 



e + 36 ( 

e v 



dx 



) = -K, 



(eqn A. 10) 



To remove the "nonlinearity" from the problem, equation (A. 10) 
is modified in the following manner: 



d 2 6 d5 

n 1 * + 3n 3 n' — - = -K (eqn A. 11) 

dx 2 dx e 
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n is defined as the approximate value of the film thickness 
from the previous iteration. In like manner, n ' is defined 
as the approximate value of the rate of change of film thick- 
ness with respect to x from the previous iteration. 

Forming the residual of equation (A. 11) yields: 

R n = n 1 * (G T d) + 3n 3 n ' (G T d) + K (eqn A. 12) 

Invoking the Galerkin criterion for the determination of the 
solution vector d, i.e. 

/G^R n dx = C i = 1,2,3,. ..n (eqn A. 13) 

yields : 

t" t 1 

nVGG ddx + 3nV/ GG ddx + K/Gdx = 0 (eqn A. 14) 

Each of the integrals in equation (A. 14) are defined in the 
following manner: 

T " n T " 

/GG ddx = Z gg d dx (eqn A. 15) 

i e 

Or on an elemental level: 

/gg T d e dx = / 

£ 

This integration results in the following 4x4 local elemental 
A matrix for any element: 



fl' 

§3 

Sa 



g l s 2 s 3 s 3 



d dx 
e 



(eqn A. 16) 
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(eqn A. 17) 



and on an elemental level 



/ ?.gg d dx = /, 



g l 

§2 

§3 

g 4 



<S 1 g 2 g 3 g 4 > 



d dx 
e 



v 

(eqn A. 18) 



This integration results in the following 4x4 local elemental 
B matrix for any element: 



[B] 



e 



-i 


z 


1 


-z 


~ T 


nr 


_ T 


nr 


- 1 
nr 


0 


z 

nr 


-£ 2 

no- 


-i 


-z 
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il 
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nr 


—j 


nr 
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nr 


z 2 

To 


- z 

nr 
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Lastly, let 



n 



/ Gdx = Z f ^ gdx 



(eqn -A . 19) 



or, on an elemental level this becomes 



I o gdx = f. 



§1 

g 2 

g 3 

g 4 



dx 



(eqn A. 20) 



This integration results in the following local elemental 
column F vector: 



l 

T 



[FJ, 



V 

n 



1 

2 

H 2 

n 



Thus, for a given element, equation (A. 12) becomes: 



[n e 4 [A] e + 3 n e 3 n e ' [B] e ] d 0 = -K e F 0 (eqn A. 21) 

As mentioned above, n and n 1 are the approximate values of 5 * 
and d<$*/dx respectively from the previous iteration. For the 
initial iteration, n is set equal to 6* ax and n* is set equal 
to 0 . 

Each elemental matrix is placed in a global system [A] 
matrix with the location in the matrix based upon the local/global 
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nodal point correspondence. For example, for element Nr. 2 
the nodal points corresponding to local nodal points 1,2,3, 
and 4 are global nodal points 3,4,5, and 6 respectively. 

Therefore, the sum of the two 4x4 elemental matrices ([A] and 
[B ] ) and multiplying constants will form a single 4x4 local 
matrix. Element (1,1) of this local matrix will be placed in 
the global A matrix location (3,3) and element (4,4) of the local 
elemental matrix will be placed in the global A matrix location 
(6,6). Ultimately, the following global system would be assembled 

[A] • d = F (eqn A. 22) 

Note that the global A matrix would be an mxm size matrix 
where m is equal to twice the number of system nodal points to 
account for the two degrees of freedom at each nodal point. 

In a similar fashion, d and F would be column vectors of size m. 

Once the system is assembled the boundary conditions are 

applied. This is done in the following manner: For boundary 

condition (a), A(l,l) is set equal to 1.0 and the remaining 

elements in the first row, i.e. A (l,i) i=2,3,4,...m are set 

equal to 0.0. Then F(l) is set equal to 6* . 6* is ini- 

m3.x max 

tially determined by a relationship developed by Leppert and 
Nimmo [Refs'. 8 and 9]. This establishes the value of d(l) as 
<5* . In a similar manner, boundary condition (b) is applied 

by setting A(2,2) equal to 1.0 and all other second row elements 
of the global A matrix are set equal to 0.0. Then, F(2) is set 
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equal to 0.0. This establishes the solution vector d(2) as 0.0. 

Note that d(2) corresponds to the slope at the first nodal point. 

One additional boundary condition is required; this is the 

value of the film thickness at x=L. This boundary condition is 

necessary to completely specify the problem. The value of the 

film thickness at the overfall may take on any value depending 

on the geometry at the overfall. In the case of this analysis, 

this value was taken as 0.25 • 6 * . This value was chosen for 

max 

the following reason: Leppert and Nimmo [Refs. 8 and 9], in 

a similar analysis for laminar film condensation on a horizontal 
surface derived an analytical solution to equation (A.l), assum- 
ing a constant surface temperature. They found the film profiles 

for 6 * with the overfall value less than 0.40 . 5 * were essen- 

max 

tially constant and thus any value of 6 * at the overfall less 

than 0.40 • 6 * would result-in the same profile. In the veri- 

fication of the finite element solution, not only was this found 

to be the case, but it was also found that the heat transfer 

rate was relatively insensitive to the shape of the film profile. 

In fact, it was found that the film thickness at the overfall 

could be increased to a value as great a 0.90 • 6 * and the 

resulting variation in heat transfer rate was only 101. This 

being the case, the value of 0.25 • 6 * was arbitrarily chosen 

max ' 

for 6 . . 

mm 

The third boundary condition was applied by setting A(m-l,i)= 
0.0 where i=l,2,3...m. Then, the global matrix element A(m-l,m-l) 
was set equal to 1.0. Finally, F(m-l) was set equal to 0.25*6* ax * 
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Once all three boundary conditions were applied, the sys- 
tem given by equation (A. 22) was solved for d. The values of 
the approximate film thickness, i.e., d(i) , i=l , 3 , 5 , . . .m- 1 are 
then compared to the values of film thickness from the previous 
iteration. If the relative difference is less than or equal 
to a specified convergence criterion (i.e., 0.0001) at all 
nodal points, convergence is met and the latest values of d 
are the solution values of 6*. 

If convergence is not met, the values of d are saved for 
the next iteration where they are used to determine n and n* 
as discussed above. This iterative process is continued until 
convergence is met or until a maximum number of iterations have 
occurred. 

B. AXIALLY FINNED CONDENSER 

The finite element solution for the film profile in an 
axially finned cylindrical condenser is very similar to that 
of a smooth cylindrical condenser. For this reason, only the 
variations in the development will be addressed. From Chapter 
II, the differential equation for mass flow rate in an axially 
finned cylindrical condenser is given by: 




5kf (T sa t - T w ) pc 
p f 2 w 2 r K fg 




) qz *1 3/ " 



(eqn A. 23) 
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This equation can be rearranged and expanded to yield: 



d 2 <5* f r*4 £*s. -v d5* f , x ±3 d6* , yi d6*-v 

(e6* 4 + 6 * tana) + (3e6* + 46* tana ^-) 



" - K ! -K 2 6* 



(eqn A. 24) 



where 1^= 3k f^ T sat ' T w^ £ 



P f 2 “ 2 r * £g 



"2 



= ? 



zcosa 



4k, CT * - T )uz* 

f_2 sat avg^ 

p, 2 u) 2 r h, cosa 
. f f § 



t 3 A 



Substituting equation (A.S) into equation (A. 24) results in 



d 2 <5 _4 5 d6 3 d6 4 d^ p 

f- e + 6 tana) + — - (3e6 - + 4T tana ^-^) 

dx 2 e e dx 6 dx e dx 



= - K, - K 7 6 
1 2 e 



(eqn A. 25) 



To remove the "nonlinearity" from the problem, equation (A. 25) 
is modified in the following manner" 



d z ~5 



dx" 



dS 

(eY 4 +Y s tana) + (3 sy 3 Y ,+ 4Y 4 Y'tana) 

dx 



= - K 1 - K 2 Y 



(eqn A. 26) 
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Here, y and y' are defined by the following relationships: 



y = 6 i + R*(^\ -6 i _ 1 ) 



(eqn A. 27) 



T 



T T 



y' = 6. + R*(S. -6 i _ 1 ) 



(eqn A. 28) 



where 6^ = approximate value of the film thickness for 
the present iteration. 

"5\ approximate value of film thickness from the 
previous iteration. 

! 

y = approximate value of the derivative of the film 
thickness for the present iteration 
75”i-l = approximate value of the derivative of the film 
thickness from the previous iteration. 

R = relaxation factor. 

These two variables, y and y', are in actuality, adjusted approx- 
imation of film thickness and derivative of film thickness re- 
spectively. This adjustment is required in order to converge to 
a solution. 

The finite element solution is now identical to that of the 
smooth condenser, that is, the residual is formed, the Galerkin 
criterion is invoked, and identical 4x4 local matrices are 
derived. Finally, the equivalent of equation (A. 21) is formed. 



t(£Y e 4 + Yg tana) [A] g + (3ey e 3 y^ + r e Y e tanC 0[B] e l • d g 



= ( - K 



le 



K 2e V le 



(eqn A. 29) 
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Notice that this equation has two forcing terms. The additional 
term (K 9 y) resulted from the nonlinearization of the problem. 

Just an in the smooth condenser film profile solution, the 
global system given by equation (A. 22) is formed, the boundary 
conditions applied and the system solved for a solution vector 
d. The iterative process is continued until convergence is met. 
With each iteration, y and y' are updated and used for the 
next iteration. When convergence is met, the latest values of 
d(i) , i=l , 3 , 5 , . . .m- 1 are the solution values of 6*. 



t 
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APPENDIX B 



USER'S MANUAL 



This appendix describes the data cards required to use 
the computer code. 

The data is divided into "blocks" for convenience. Each 
page of this user's guide is a separate block. For each block, 
a general description, the required format, and appropriate 
comments are provided. 

It is imperative that input data be consistent or errors 
will result. For example, if a smooth geometry is being analyzed, 
the finite element parameters must also result in a smooth model. 
In addition, all data fields must be filled with an input value, 
even if that value is not needed for the analysis. For example, 
in a smooth analysis, no fin half angle is required for the cal- 
culations; however, a value of the correct format must be pro- 
vided in the fin half angle field or an INPUT/OUTPUT error will 
result . 
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DATA BLOCK A 

DESCRIPTION: FINITE ELEMENT PARAMETERS 

FORMAT : 815 

12345678 



NDIV NCMREC NCMTRI NRWFIN NRWTRF NCMTRF NCOL NPRNT 



FIELD CONTENTS 

1 NDIV Number of axial increments. Must be less than 

or equal to 50. 

2 NCMREC- Number of columns of finite elements in the rec- 

tangular section of fin. May be equal to 0 if 
fin is triangular only. 

3 NCMTRI-Number of columns of finite elements in tri- 

angular section of fin. May be equal to 0 if 
only rectangular fin. 

4 NRWFIN- Number of rows of finite elements in the fin 

section of model. May be equal to 0 if a smooth 
condenser . 

Note: If triangular or trapezoidal fin, NRWFIN 

must equal NCMTRI. 

5 NRWTRF- Number of rows of finite elements in wall section 

of the model. 

6 NCMTRF-Number of columns of finite elements in the 

trough section of finned model. Set equal to 
0 if a smooth condenser. 

7 NCOL — Total number of columns of finite elements. 

Must be equal to NCMREC+NCMTRI+NCMTRF for a 
finned model. 

8 NPRNT- -Print control number. 

If equal to 1 correspondence table and major 
elements of finite element model will be printed 
out . 

If equal to 0, output will be suppressed. 
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DATA BLOCK B 



DESCRIPTION: FLUID, FIN MATERIAL SELECTION PARAMETERS 



FORMAT : 


215 


1 


2 3 4 5 6 7 8 


IFLUID 


IFIN 



FIELD 


CONTENTS 


1 


IFLUID-If equal to 0, working fluid is water. 

If equal to 1, working fluid is freon. 


2 


IFIN If equal to 0, condenser wall material is 

copper . 

If equal to 1, condenser wall material is 
stainless steel. 
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DATA BLOCK C 



DESCRIPTION: CONDENSER GEOMETRY 

FORMAT: 6G10.5 

12345678 



CLI REASEI THICKI BFINI CANGL FNWTHI 



FIELD CONTENTS 

1 CLI ---- Condenser length (inches). 

2 RBASEI - Inside radius to wall of condenser at condenser 

end (inches ) . 

3 THICKI-Wall thickness (inches). 

4 BFINI--Fin height (inches). Must be set equal to 

0.0 if smooth condenser. 

5 CANGL-- Condenser half angle for tapered condenser 

(degrees). Must be set equal to 0.0 for 
cylindrical condenser. 

6 FNWTHI-Width of rectangular portion of trapezoidal 

or rectangular fin (inches). If triangular 
fin, set equal to 0.0. 
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DATA BLOCK D 



DESCRIPTION: INTERNAL FIN GEOMETRY 



FORMAT : 


2G10 . 5 


1 


2 3 4 5 6 7 8 


FANGL 


ETOEO 



FIELD 


CONTENTS 


1 


FANGL- -Fin half angle (degrees). Set equal to 0.0 
for smooth condenser or rectangular fin. 


2 


ETOEO--Ratio of trough width to fin base width. 
Determines spacing between fins. 
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DATA BLOCK E 



DESCRIPTION: CONVERGENCE CRITERIA 



FORMAT: 


2G10 . 5 


1 


2 3 4 5 6 7 8 


CRIT 


CRITDL 



FIELD 


CONTENTS 


1 


CRIT Temperature convergence criterion. Used to 

determine solution of two dimensional steady 
state heat conduction problem. 


2 


CRITDL-Mass flow convergence criterion. Used only 
in cylindrical condenser analysis for mass 
flow convergence test. 
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DATA BLOCK F 



DESCRIPTION: OPERATING PARAMETERS 

FORMAT: 5G10.5 

12345678 



RPM HINF TINTL TSAT TINF 



FIELD CONTENTS 

1 RPM Rotational speed (revolutions per minute). 

2 HINF External heat transfer coefficient 

(Btu/hr-ft -deg F) . 

3 TINTL- - Initial temperature estimate (degrees F) . 

4 TSAT Saturation temperature (degrees F) . 

5 TINF Ambient temperature (degrees F) . 
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DATA BLOCK G 



DESCRIPTION: OUTPUT PRINT CONTROL 



FORMAT : 


415 


1 


2 3 4 5 6 7 8 


IUNITS 


NFLAG1 NFLAG2 NFLAG3 



FIELD 


CONTENTS 


1 


IUNITS -Output units control number. 

If IUNITS = 0, calculated results will be pro- 
vided in English units. 

If IUNITS = 1, input parameters will be repeated 
in SI units and calculated results will be pro- 
vided in SI units. 

If IUNITS = 2, input parameters will be repeated 
in SI units and output results will be provided 
in both English and SI units. 


2 


NFLAGl-The first axial increment at which the parameters 
listed under remarks will be provided as output. 


3 


NFLAG2-The final increment at which the parameters 
listed under remarks will be provided. 


4 


NFLAG3-The step change in increments between NFLAG1. 


REMARKS 





1) No matter what value of IUNITS is used, input parameters 
will always first appear in English units. 

2) For increments indicated by NFLAG values, the following 

parameters will appear: a) x- coordinate , y-coordinate and 

temperature at each nodal point, b) length of element, heat 
transfer coefficient and heat rate per unit length for each 
convective boundary element. 

3) As a minimum, the values of 1, 1, 1 must be provided as 
input for NFLAG1, NFLAG2 , and NFLAG3 respectively. 
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DATA BLOCK H 



DESCRIPTION: TAPERED CONDENSER SOLUTION METHOD 

FORMAT: 115 

12345678 



NSOLVE 



I 



FIELD CONTENTS 

1 NSOLVE-Tapered solution control number. 

For tapered, axially finned condenser, set 
NSOLVE=l . 

For tapered smooth condenser, NSOLVE must be 
set to one of the following three values: 

Set NSOLVE = 2 if solution of film thickness is 
to be based on Ballback's [Ref. 1] 
equation. 

Set NSOLVE = 3 if solution of film thickness is 
to be based on Daniels and Al-Jumaily [Ref. 13] 
equation, neglecting drag terms. 

Set NSOLVE = 4 if solution of film thickness is 
to be based on Daniel's and Al-Jumaily [Ref. 

13] equation, with drag effects included. 



REMARKS 

1) NSOLVE is only used in tapered condenser analysis. 



115 



DATA BLOCK I 



DESCRIPTION: CYLINDRICAL CONDENSER ANALYSIS PARAMETERS 



FORMAT : 


315 ,2GI5. 10 


1 


2 3 4 5 6 7 8 


NONCE 


ITERMX ITRPRT RELAX DELMAX 



FIELD 


CONTENTS 


1 


NONCE- - Single iteration parameter. 

If NONCE = 1, only one iteration will be 
permitted . 

If NONCE = 0, iterations will continue until 
convergence or maximum number of iterations 
is reached. 


2 


ITERMX-Maximum number of iterations permitted in analysis. 


3 


ITRPRT- Iteration print control parameter. 

If ITRPRT = 1, mass flow convergence test results 
will be provided for each iteration. 

If ITRPRT = 0, mass flow convergence test results 
will only be provided on final iteration. 


4 


RELAX- -Relaxation variable used in finite element solution 
of cylindrical finned film profile. 


5 


DELMAX- Initial estimate of maximum film thickness used in 
solution of cylindrical finned film profile. 


REMARKS 





1) Above parameters are only used in cylindrical condenser 
analysis . 

2) Recommended value of RELAX is 0.80. It is sometimes necessary 
to adjust this value plus or minus 0.05 to reach film profile 
convergence at small film thickness values. 

3) Input value of DELMAX is only used in cylindrical finned 
analysis. For the cylindrical smooth condenser, DELMAX is 
internally generated. 
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APPENDIX C 



SOURCE CODE LISTING 
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INPUT SCLUT10N METHOD VARIAdLE FOR TAPERED SMOOTH CASE 
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OUT PL T MAJOR TEMPERATURES FOR EACH INCREMENT 
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